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ABSTRACT 

We present the analysis of high-resolution optical spectroscopic observations of the zero-age main-sequence O 
star Herschel 36 spanning six years. This star is definitely a multiple system, with at least three components 
detected in its spectrum. Based on our radial-velocity (RV) study, we propose a picture of a close massive 
binary and a more distant companion, most probably in wide orbit about each other The orbital solution for 
the binary, whose components we identify as 09 V and BO. 5 V, is characterized by a period of 1 .5415 ± 0.0006 
days. With a spectral type 07.5 V, the third body is the most luminous component of the system and also 
presents RV variations with a period close to 498 days. Some possible hypotheses to explain the variability are 
briefly addressed and further observations are suggested. 

Subject headings: binaries: close — stars: early-type — stars: individual (Herschel 36) 



1. INTRODUCTION 

The 0-type star Herschel 36 (Her 36) is located in the 
M8 Nebula, a high-mass star-forming region at a distance of 
1.3 kpc (Arias et al. 2006). This 07.5 V(n) star (Walborn 
1982) is responsible for the ionization of the optically bright- 
est part of M8, known as the Hourglass Nebula. A very young 
infrared cluster is found around Her 36 (Arias et al. 2006). 
The ultracompact Hll region G5. 97-1. 17, a candidate "pro- 
plyd" similar to those seen in the Orion Nebula (Stecklum 
et al. 1998), is located 2'.'7 distant from Her 36. Recently, 
high-resolution infrared imaging of the region also revealed 
the existence of a compact source at 0'.'25 southeast of the 
star, which harbors an early B-type embedded star (Goto et 
al. 2006). Her 36 and its close neighbors are probably part of 
a Trapezium-like stellar system. 

It seems clear that most, or even all, massive stars form 
in clusters and OB associations (de Wit et al. 2005). After 
a few hundred thousand years, they become optically visi- 
ble. Only few, if any, of the massive stars are found as sin- 
gle objects and, on average, massive stars have more than one 
companion. Ignoring the multiplicity of massive stars can in- 
troduce serious problems when comparing observations and 
stellar evolution theory. 

Although it is often stated that zero-age main-sequence 
(ZAMS) O stars should not be and are not observed optically. 
Her 36 appears as a possible bona fide ZAMS O star (Wal- 
born 2007), based on the appearance of its optical spectrum, 
the association with a dense, dusty nebular knot and its ob- 
served subluminosity (Arias et al. 2006). Because of all these 
aspects. Her 36 and the Hourglass have been many times com- 
pared to 9^ Ori C and the Orion Trapezium. Her 36 has also 
been extensively studied because it shows an "anomalous" ex- 
tinction curve with Ry = 5.36 (Arias et al. 2006), which rep- 
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resents one of the highest Ry values known (cf. Fitzpatrick 
& Massa 2009). 

A long-term spectroscopic monitoring of Galactic O and 
Wolf-Rayet of the Nitrogen sequence (WN) stars is being car- 
ried out (Gamen et al. 2008), with the main aim of estab- 
lishing the multiplicity status of these objects. As part of this 
observational campaign. Her 36 has been systematically ob- 
served. Here we present the results of the analysis of the data, 
which demonstrate conclusively that Her 36 is a complex mul- 
tiple system, with outstanding characteristics. 

2. OBSERVATIONS AND DATA REDUCTION 

Table[T]presents the journal of the observations used in this 
work. Most of them came from the monitoring campaign of 
Galactic O and WN stars and were obtained between 2005 
April and 2009 May at the observatories of Las Campanas 
(LCO) and La Silla, in Chile. A total of 33 high-resolution 
spectra were obtained at LCO; 32 of those were taken with 
the echelle spectrograph at the 2.5 m Du Pont telescope and 
the other with the MIKE at the 6.5 m Magellan II (Clay) tele- 
scope. Twenty-six more spectra were acquired with FEROS at 
ESO La Silla. Four additional spectra were taken at the Com- 
plejo Astronomico El Leoncito (CASLEO), Argentina, in ob- 
serving runs prior to 2005, using the REOSC SEL Cassegrain 
spectrograph in cross-dispersion mode. Technical details of 
the observations and data reduction are presented by R. H. 
Barba et al. (2010, in preparation). 

3. THE OPTICAL SPECTRUM 

Besides the Balmer series, the spectrum of Her 36 shows 
many absorption lines characteristic of hot massive stars, such 
as those of He I, He II, C III and IV, O III and Si IV. Some dif- 
fuse intersteflar bands at 4430, 4500, 5780, and 5797 A, as 
well as a number of interstellar absorption lines due to Ca II 
and Na I, are also evident. Narrow strong nebular emission 
lines, associated with the Hourglass Nebula, are observed. 
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some of them superimposed on the Balmer and He I absorp- 
tions. 

The spectrum of Her 36 was monitored during two to six 
consecutive nights in each of the observing runs. The first 
result that came from the inspection of our high-resolution 
data is that the spectrum of this star is highly variable. 

The spectral lines of Her 36 show large variations in both 
line-profile shape and radial velocity (RV), on a timescale 
of hours. Even if all the absorption lines appear apprecia- 
bly broadened, in good agreement with the qualifier (n) in the 
spectral classification by Walborn (1982), the profiles of the 
most prominent ones alternate between more or less normal 
and unusually broad. When observed with enough resolution, 
the lines are clearly not single: at different times and in differ- 
ent lines, one, two, or three components are detected, showing 
that this star is actually a multiple system. 

We present the variations in the He II A4686 line in the left 
panel of Figure [T] This plot includes the data from five con- 
secutive nights. The morphology observed varies from a very 
deep single profile to a well-defined double line. At first sight, 
however, the line-profile morphology does not seem to be 
strictly repeatable. Large RV variations accompany all these 
changes. We will refer to the main and secondary absorptions 
as "component A" and "component Bl", respectively. 

An even more complex behavior is observed in the He I ab- 
sorption lines. The right panel of Figure [T] shows the spectral 
variations in the He I A5875 line. Now three absorption com- 
ponents (labeled as A, Bl, and B2) can be clearly identified. 
Whereas components A and Bl correlate with those observed 
in the He II A4686 profile, the component B2 is undetected 
in the He II lines, indicating that the latter star is cooler than 
its companions. The RV of the main component A changes 
rather slowly, but the variations observed for components Bl 
and B2 are remarkably large. Note that the positions of the 
components B 1 and B2 practically invert from the forth to the 
fifth night. 

Apart from the He II A4686 and He I A5875 lines, several 
other photospheric features, such as He I A4471, He II A5411, 
Cm A4650, and Civ A5812, were also measured. We found 
correlated variations in profile shape and RV in most of these 
lines. 

Finally, we considered the least blended spectra in order to 
obtain a spectral type for each of the three stellar components 
detected in the spectrum of Her 36. The spectral classifica- 
tion was performed by comparison of the spectrum with the 
OB star adas of Walborn & Fitzpatrick (1990). We adopt the 
spectral types 07.5 V, 09 V, and B0.5V for components A, 
Bl, and B2, respectively. 

4. A PICTURE FOR THE MULTIPLE SYSTEM HER 36 

The fact that all measured photospheric lines show similar 
profile and RV variations suggests that the RV variations de- 
tected in the spectrum of Her 36 are actually caused by the 
orbital motion of, at least, three stellar components. We pro- 
pose that the system is composed of a close O-type binary and 
a third more distant companion, also of spectral type O. The 
binary and the third body are very likely in wide orbit about 
each other. 

4. 1 . The Close Pair Bl +B2 

The behavior observed in the He I A5875 line suggests that 
the components B 1 and B2 are constituents of a double-lined 
spectroscopic binary. When we refer to this binary system we 
will simply use the term "system B". 



As mentioned before, the component B2 is not detected in 
He II. Moreover, only in the He I A5875 line, the profile fea- 
tures appear separated enough to obtain a reliable measure- 
ment of the RV of each component. Thus we considered only 
this line to determine the orbital solution for the stellar system 
B. The RV measurements used in the determination of the or- 
bital solution are presented in Table [T] (Columns 6 and 7). We 
determined the RVs of the different components of the He I 
A5875 line by fitting them with adequate sums of Gaussian 
profiles. To do this, we considered four Gaussian functions: 
one for each of the stellar components A, Bl, and B2, as well 
as an additional one to take into account the strong nebular 
emission. First, we used the least blended spectra to deter- 
mine the intensities and widths of the three stellar blends and 
then repeated the disentangling by keeping these parameters 
fixed. Thus the deconvolution of Gaussian profiles yielded the 
central wavelengths of the four components (three stellar and 
one nebular) of the He I A5875 line as adjusted parameters. 
In spite of the complexity of this blend, the fits were highly 
satisfactory in most of the cases. Being conservative, we es- 
timate that the measurement errors in the RVs of the stellar 
components are no larger than 10 kms^^. 

To derive the period of the orbital motion, we applied a va- 
riety of methods to the RV measurements corresponding to 
the components Bl and B2. These routines include the Lafler 
& Kinman (1965) method and subsequent modifications of 
it (Man-aco & Muzzio 1980, MM80), as well as information 
entropy minimization (Cincotta et al. 1995, CMN95). We 
considered three RV data sets: (1) the measurements of the 
primary component, Bl, alone; (2) the measurements of the 
secondary component, B2, alone; and (3) the measurements 
of both components together. Based on the resulting peri- 
odograms, we adopted P — 1.54 days as the best estimate for 
the period. 

To derive the orbital elements of the system B, we used 
a modified version of the code originally written by Bertiau 
& Grobben (1969), adopting the value Pq = 1.54 days as 
starting estimate of the period and considering the same data 
sets as previously described. After some trials with general 
eccentric fits, a circular orbit was assumed for the system. The 
best orbital solution is presented in Tableland illustrated by 
the RV curves in the upper panel of Figure |2] 

4.2. The Most Luminous Star A 

The star A is the most luminous component of the system 
and dominates the optical spectrum of Her 36. Although not 
so evident, it also presents RV variations. Careful inspection 
of the data suggests the presence of a long-term (few hundred 
days) variability. 

To study the spectral variations, we used the RV measure- 
ments of the He II A4686 line. This line appears as one of the 
most useful features for investigating possible periodicities 
because of being less contaminated by other features or the 
background nebular emission, and, mainly, because of show- 
ing only two stellar components: A and Bl. As in the case 
of the He I A5875 line, we determined the RVs of the differ- 
ent components by performing multiple Gaussian fits. These 
measurements are presented in Columns 3 and 4 of Table [T] 
Due to the strong overlap of the line components, our initial 
disentangling attempts yielded systematic shifts of the central 
position of the component A, which at first we erroneously in- 
terpreted as short-term periodic RV variations. However, this 
effect was fully connected when we repeated the disentangling 
by keeping the intensities and widths of the Gaussian pro- 
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Fig. 1. — Variations in tlie Hell A4686 (left panel) and He I A5875 (right panel) line profiles, observed in LCO/Du Pont spectrograms obtained in five 
consecutive nights. The spectra are labeled with the corresponding Juhan Date (JD) in the format HJD-2,450,000. All the spectra are normalized to unity and 
have equal scales. The separation between continuum levels is 10% of the normalized flux. A, Bl, and B2 refer to the components described in the text. The 
strong emission line superimposed on the He I absorption is associated with the Hourglass Nebula. 



TABLE I 

Radial-velocity measurements of the stellar components observed in the spectra of Her 36 



No. 


HJD 


He II 4686 




He I 5876 




Telescope 


Instrument 




-2450000 


A 


Bl 


A 


Bl 


B2 






1 


2471.6791 


-10.2 


126.0 








Casleo-2.15 


REOSC 


2 


2835.6685 


-58.8 


92.9 








Casleo-2.15 


REOSC 


3 


2836.7348 


-72.5 


218.4 


-65.0 


256.8 


-229.5 


Casleo-2.15 


REOSC 


4 


2837.6846 


-54.5 


-146.9 








Casleo-2.15 


REOSC 


5 


3490.8852 


26.4 


-128.7 








LCO-DP 


Echelle 


6 


3491.9145 


13.6 


205.3 


16.6 


213.8 


-336.1 


LCO-DP 


Echelle 


7 


3491.9199 


19.3 


210.5 


18.6 


231.4 


-338.8 


LCO-DP 


Echelle 


8 


3772.8789 


-47.0 


50.3 








LCO-DP 


Echelle 


9 


3873.9131 


-50.9 


91.2 








LCO-DP 


Echelle 


10 


3874.9028 


-44.9 


-175.5 


-46.9 


-132.5 


347.3 


LCO-DP 


Echelle 



Note. — Rest wavelengths used in the determination of the RVs are 4685.65 A and 5875.65 A, respectively. 
(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 



files fixed, in order that the wavelength positions of A and B 1 
were the only parameters adjusted in the fit. As a result, we 
obtained very good quality fits and reliable RVs, with typical 
measurement errors no larger than 10 km s^^. 

To estimate the possible period of the variability of A, we 
applied the same period-finding routines as in Sec. 14.11 i.e., 
the methods from MM80 and CMN95, to the correspond- 
ing RV measurements, and we also performed the Lomb- 
Scargle (LS) periodogram search (Scargle 1982). A num- 
ber of possible values of the period can be identified from 
the periodograms peaks, the most relevant corresponding to 
P = 508 ± 15 days. An additional value of P 206 days 
is also strong but less significant, and cannot be discarded a 
priori. 



Considering the values from the periodogram peaks as 
starting estimates of the period, we explored the possibility of 
deriving an orbital solution for the motion of the component 
A. For this purpose, we applied the same code as in Sec. 14. II 
to the RVs of the component A of the Hen A4686 line. 
The resulting solution is presented in the second column 
of Table |3] and illustrated by the RV curve in the middle 
panel of Figure |2] An independent solution based on the RV 
measurements of the component A of the He I A5876 line 
(Column 5 of Table [Til was also calculated, and it is presented 
in the third column of Table |3] It can be verified that, within 
the errors, both solutions represent essentially the same orbit. 
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Fig. 2. — Upper panel: RV curves coiTesponding to the close stellar pair 
B1+B2. The filled and open circles indicate the observed RVs of the primary 
(Bl) and secondary (B2), respectively. Zero phase is a time of maximum 
velocity for the primary Bl. Middle panel: the RV curve derived from the 
measurements of the strongest component (component A) of the He 11 A4686 
line. Lower panel: RV curves corresponding to a joint orbital solution ob- 
tained from the He I A5876 fine, considering the RVs of the component A for 
the primary (filled squares) and the systemic velocities of the system B for 
the secondary (open squares). The orbital parameters that characterize these 
solutions are fisted in Table[3] Zero phase is a time of maximum velocity for 
the component A. 



TABLE 2 

Orbital and physical parameters of the system B 



Ps (days) 

Tvmax (HJD-2,450,000) 

Vo (kms-1) 

Ki (kms-l) 

K2 (kms-l) 

ai sin is (km) 

a2 sin is (km) 

Ml sin^ ig (Mq) 

M2 sin^ is (M©) 

Q (M2/M1) 

rms (km s~^) 



1.5415 ±0.0006 
4248.8053 ± 0.0061 
6.3 ±3.6 
201.4 ±5.8 
327.3 ±5.8 
(42.7 ± 1.2) X 10^ 
(69.4 ± 1.2) X 10^ 
14.6 ± 1.3 
9.0 ±1.3 
0.61 ± 0.03 
31.0 



Note. — Parameters derived from the RVs of the He I A5875 line and 
assuming a circular orbit. Tvmax refers to the maximum velocity of the 
primary star B 1 . 



5. SUMMARY AND DISCUSSION 

The main finding of this study is the multiple nature of the 
ZAMS object Her 36. The absorption lines in the spectrum of 
this star are clearly not single. We arrived at the conclusion 
that Her 36 presents at least three stellar components that we 
labeled as A, Bl, and B2, and classified as 07.5 V, 09 V, and 
BO. 5 V, respectively. 

Arias et al. (2006) determined the distance and reddening 
to several early-type stars in the Hourglass region, including 
Her 36. As mentioned in Sec.[T] Her 36 that was considered 
as a single 07.5 V star in their calculus fell below the calibra- 
tion of its luminosity class. The here demonstrated multiplic- 
ity of Her 36 makes the stars that compose this triple system 
even much more subluminous individually. Thus, if we adopt 
the absolute magnitudes for ZAMS OB stars from Hanson 
et al. (1997), the combined My value obtained is similar to 
that quoted by Arias et al. (2006). This fact adds additional 
support to the hypothesis that the stellar components of this 
system may be on the ZAMS. 

We measured the most prominent features in the spectrum 
of Her 36, finding correlated variations in profile shape and 
RV. Based on our RV study, we proposed that Her 36 is com- 
posed of a close massive binary, the stellar pair B1-I-B2 (sys- 
tem B), and a companion, the component A, which is the most 
luminous star of the system and dominates the optical spec- 
trum. 

Considering the RV measurements of the He I A5875 line, 
we derived a set of orbital elements for the system B. The 
best solution is characterized by a circular orbit with a pe- 
riod Ps — 1.5415 ± 0.0006 days, and semi-amplitudes Ki = 
201.4 kms^i and K2 = 327.3 kms"^ We can use the ob- 
tained minimum masses M sin^ i to estimate the inclination 
of the system B. The values of 14.6 and 9.0 Mq, derived 
for the primary and the secondary, respectively, are compared 
with the typical masses of 09 V and BO. 5 V stars. According 
to the values for detached eclipsing systems from Gies (2003) 
(19.0 Mq for 09 V and 1 1 .7 Mq for B0.3 V), the system must 
have a large inclination (larger than 67°), which means that 
it could possibly undergo eclipses. Unfortunately, no photo- 
metric study of Her 36 exists at present. Future photometric 
monitoring will certainly be valuable. 

The component A also shows RV variations. Independent 
RV measurements of the He II A4686 and the He I A5876 lines 
yield essentially the same relatively eccentric (e = 0.3) or- 
bit, characterized by a period close to 498 days and a semi- 
amplitude of the order of 65 km s ^ ^ . 

We feel that the component A and the close binary system B 
are in wide orbit about each other. The long-term periodicity 
observed in the RVs of the component A must then be related 
to the mutual orbit of both objects. If that is the case, and we 
are indeed in the presence of a triple (hierarchy 2-1-1) system, 
then we should expect the systemic velocity of the binary B to 
change with a period Pl ~ 498 days. To prove this hypothe- 
sis, we computed a joint orbital solution considering the RVs 
of the component A of the He I A5876 line for the primary, 
and the systemic velocities of the system B, vb, for the sec- 
ondary. The Vb values were calculated according to the basic 
expression Vg = {MiVi + M2V2)/ {Mi + A/2), where Mi and 
M2 are the minimum masses from Table |2l and vi and V2 are 
the RVs of the components Bl and B2 of the He I A5876 line, 
respectively. The resulting solution is characterized by the pa- 
rameters in the last column of Table [3] and the RV curves in 
the lower panel of Figure |2l Although with somewhat larger 
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dispersion, these curves represent practically the same orbit 
as those derived from the RV data of A alone, as can be easily 
verified by confronting the orbital parameters listed in the dif- 
ferent columns of Table|3] We note here that the joint solution 
was obtained in a completely independent manner, supporting 
the idea of the mutual orbital motion between A and B. 

Comparing the projected mass of the system B from Ta- 
ble [3] with the projected total mass of the close pair B1H-B2 
from Table|2] it can be shown that sini^ = 1.03 sinis- This 
means that the inclinations of the two orbital planes satisfy 
either is ~ m or is ~ 180° — iA- In the first case, the inner 
and the outer orbits of the triple system will be nearly copla- 
nar, whereas in the last case the relative inclination between 
the orbital planes, 0, will He between 180° - 2 < </> < 180° 
(Fekel 1981). The relative orientation of the orbits in a mul- 
tiple system is of particular interest because of its direct rela- 
tionship with the conditions at the epoch of formation. Thus, 
the study of a triple system as young as Her 36 can provide 
important clues to the formation processes of massive stellar 
systems. 

On the other hand, comparing the minimum mass for the 
primary A (19.2 Mq) to typical masses of 07.5 V stars from 
Gies (2003), we estimate that iA is most probably larger than 
70°. The component A has the largest individual mass of the 
three stars. Its mass is however fairly comparable with the 
total mass of the inner binary B. According to Table [3] the 
minimum mass of the triple system is 19.2+26.0 = 45.2 Mq. 



Assuming a minimum inclination ia — 70°, we conclude that 
the total mass of Her 36 is probably no larger than 54 Mq. 

The triple ZAMS system Her 36 is really peculiar, and it 
undoubtedly deserves additional study. Clearly, further obser- 
vations are required before a complete characterization of the 
orbits can be given and a definitive model for the system can 
be derived. Spectroscopic observations with higher signal-to- 
noise ratio, as well as an exhaustive program of photometric 
monitoring, seem essential. 
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TABLE 3 

Orbital and physical parameters of the multiple system Her 36 



Parameter HellA4686 HelA5876 HelA5876 
A A A+B 

Pl (days) 498.3 ± 3.7 496.3 ± 2.2 493.3 ± 3.2 

e 0.29 ±0.04 0.30 ± 0.03 0.25 ± 0.03 

a)(deg) 126.0 ±6.3 126.9 ± 3.4 124.5 ± 5.4 

To (HJD-2,450,000) 4279.5 ± 7.5 4278.6 ± 3.9 4275.6 ± 6.5 

Tvmax (HJD-2,450,000) 4147.5 ± 7.5 4146.7 ± 3.9 4141.5 ± 6.5 

\/o(kms-i) 5.3 ±1.8 4.4 ± 0.9 7.4 ± 0.8 

_ft:A(kms-i) 65.2 ±3.1 61.2 ± 2.2 57.2 ±2.4 

iiTs (km.s-1) - - 12.1 ±2.4 

aAsinjA(km) (42.7 ± 3.0) X lO'' (39.8 ± 1.9) X lO'^ (37.5 ± 2.1) x lO'' 

as sin (km) - - (27.6 ± 2.0) X lO'' 

/(M) (Mq) 12.5 ±2.4 10.2 ±1.5 

Ma sin3jA(M0) - - 19.2 ±6.1 

Mgsin^ iA(^Q) - - 26.0 ±6.3 

Q{Mb/Ma) - - 1.36 ±0.13 

rms(kms-i) 10.3 3.6 7.0 



Note. — Columns 2 and 3 present the orbital parameters derived from the RVs of the component A alone, considering the He II A4686 and the He I A5876 
lines, respectively. Column 4 presents a joint orbital solution obtained from the He I A5876 Une, considering the RVs of the component A for the primary and the 
systemic velocities of the system B for the secondary. 



